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By considering a flavour expansion about the St/ (3)-flavour symmetric point, we investigate how 
jy-^ flavour-blindness constrains octet baryon matrix elements after SI/ (3) is broken by the mass dif- 

ference between the strange and light quarks. We find the expansions to be highly constrained 
along a mass trajectory where the singlet quark mass is held constant, which proves beneficial for 
extrapolations of 2+1 flavour lattice data to the physical point. We investigate these effects nu- 
merically via a lattice calculation of the flavour-conserving and flavour-changing matrix elements 
of the vector and axial operators between octet baryon states. 
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1. Introduction 

Understanding the pattern of flavour symmetry breaking and mixing, and the origin of CP 
violation, remains one of the outstanding problems in particle physics. In [1, 2] we have outlined a 
program to systematically investigate the pattern of flavour symmetry breaking. The program has 
been successfully applied to meson and baryon masses involving up, down and strange quarks. In 
these talks we will extend the investigations to include matrix elements. 

The QCD interaction is flavour-blind. Neglecting electromagnetic and weak interactions, 
the only difference between flavours comes from the mass matrix. We investigate how flavour- 
blindness constrains matrix elements after flavour SU (3) is broken by the mass difference between 
the strange and light quarks, to help us extrapolate 2 + 1 flavour lattice data to the physical point. 

We have our best theoretical understanding when all 3 quark flavours have the same masses 
(because we can use the full power of flavour SU (3)); nature presents us with just one instance of 
the theory, with m s /mi « 25. On the lattice we can choose our quark masses, so we can investigate 
fictional universes where m s /mi / 25, and so gain a clearer understanding of flavour symmetry 
breaking. 

We have previously used a symmetry analysis, similar in spirit to that used by Gell-Mann 
and Okubo [3, 4] in the earliest days of the quark picture, to find formulae for the quark mass 
dependence of hadron masses [1, 2]. We now extend this analysis to hadron matrix elements. In 
the first part of these proceedings we discuss the group theory, in the second part we compare our 
expectations with lattice data. We shall then discuss briefly the further steps needed to compute 
form factors relevant to the determination of the CKM matrix element \V US \- 

In this work we concentrate on the 2 + 1 case, in which symmetry breaking is due to mass 
differences between the strange and light quarks; but our methods are also applicable to isospin 
breaking effects coming from a non-zero — m u , e.g. [5]. 

2. St/ (3) breaking 

How severely does the strange quark mass break SU (3) symmetry? In this approach it is not 
the strange-light ratio, m s jm\ ~ 25, which matters. A more natural way to judge the severity of 
symmetry breaking is to compare (m s — mi) with a typical hadronic mass. Since (m s — m{) <C Mg 
we can hope for an expansion with very good convergence. 

We can see how well this works in practice by looking, for example, at the physical masses 
of the decuplet baryons. We can construct mass combinations which first appear at different or- 
ders of the symmetry breaking, starting with quantities that would be non-zero even with perfect 
SU (3), and working up to quantities which first appear at the third order in the symmetry breaking 
parameter 



4M A + 3M Z * + 2M Z * + Mq, 
-2M A +M s *+M n 
4M A - 5M E » - 2M Z * + 3Afo 
-M A + 3M Z . - 3M S * +M& 



+0.742 GeV 



+ 13.82 GeV 



-0.044 GeV 



-0.006 GeV 



27-plet 
64-plet 



octet 



singlet 




(2.1) 
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where we have used the isospin-averaged experimental masses for each decuplet baryon. Clearly 
we see a strong hierarchy in values. Each additional factor of (m s — mi) reduces the value by about 
an order of magnitude, the final 0((m s — mi) 3 ) quantity is about 2000 times smaller than the leading 
quantity. An expansion that yields a factor of 10 for each order is very good compared with most 
approaches available for QCD. 

To investigate flavour symmetry breaking systematically, we need to vary the amount of sym- 
metry beaking we have, while keeping all the flavour singlet terms in the action constant. We 
therefore follow a strategy in which the average quark mass m = (m u + m d + m s )/3 is kept con- 
stant, while the mass splitting is increased. Our notation for the quark masses and their splittings 
is 

m = -(m u +md + m s ) fixed 
8m u = m u — m 

8m d =m d -m (2.2) 
8m s = m s —m 

mi = -{m u +m d ) 
8mi = mi —m. 

From these definitions we have the identity 8m u + 8m d + 8m s = 0. In this notation the quark mass 
matrix is 

Jt = 



m u 











m d 











m s 









( : 


1 




Vo 













°\ 






-1 





| + \8m s 


Vo 





o) 






ml 1 I +\{8m u -8m d ) I -1 0]+|5mJ -1 I . (2.3) 
\0 1/ \0 0/ \ 2/ 

^# has a flavour singlet part (proportional to /) and a flavour octet part, proportional to A3, X%. It 
is important to note that there are no terms in the QCD Lagrangian which are in representations 
higher than the octet. The only way to give a value to a quantity in a higher SU (3) representation 
is to have multiple powers of the flavour-breaking term, i.e. multiple powers of 8m q . 

Therefore we adopt the following strategy. We classify physical quantities by their representa- 
tion of SU (3) and its sub-group SU (2), and classify quark mass polynomials in the same way. The 
Taylor expansion of a quantity of known symmetry can only involve polynomials of the matching 
symmetry. This strongly constrains the Taylor expansion of physical quantities about a symmetric 
point with all three quark masses equal. 

In this work we are investigating non-singlet matrix elements (e.g. vector and axial-vector 
currents for weak decays) acting between octet baryons. This octet is illustrated in Fig. 1. The 
matrix elements have the form (B' | \B ) with i, j, k G {1, ■ ■ • , 8} because the hadrons and operators 
are both in flavour octets. Our choice of indices is set out in Table 1. We use the corresponding 
meson name to refer to the flavour of bilinear quark operators, for example 

(p\n°\p) = ( p \±=(uyu-dyd)\p) = (B 2 \0 4 \B 2 ), (2.4) 
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Figure 1: The (lowest) octet baryon states. 



Index 


Baryon (5) Meson (M) Operator (0) 


1 

2 
3 
4 
5 
6 
7 
8 


7i #° dys 
p K + uys 

1° 71° ^(uyu-dyd) 
A ^ 7S ("^ M + ~ 

S~ i<r syu 



Table 1: Our numbering of the octet states, used internally. Whenever possible, final results will be presented 
in ways that are independent of the index choice. We use the convention that operator number i has the same 
effect as absorbing a meson with the index i. 



where 7 is a generic gamma matrix. Since we need three indices, i, j, k, to specify a matrix element, 
we now need to classify 8x8x8 tensors under SU (3), in just the same way as we needed to give 
the classification of 8 x 8 and 10 x 10 matrices for baryon masses. 

To find the allowed mass-dependence of octet matrix elements of octet hadrons we need the 
SU (3) decomposition of 8 8 8. Using the intermediate result 

8®8 = 1©8©8©10©I0©27, (2.5) 

we find 

80808=l©lffi8ffi8ffi8ffi8ffi8ffi8ffi8©8 
©27 ©27 ©27 ©27 ©27 ©27 ©64 
ffilOffilOffilOffilOffilOffilOffiTOffiTO 
035 35 35 35. (2.6) 

The allowed quark mass Taylor expansion for a hadronic matrix element must follow the 
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schematic pattern 

(H^O^H*) = £(singlet mass polynomial) x (singlet tensor) l] 

+ ^(octet mass polynomial) x (octet tensor)' 7 * (2.7) 
+ ^(27-plet mass polynomial) x (27-plet tensor)' 7 * 
+ ^(64-plet mass polynomial) x (64-plet tensor)' 7 * 
+ ■■■ . 

The tensors in this equation are three-dimensional arrays of integers and square-roots of integers, 
objects somewhat analogous to three-dimensional Gell-Mann matrices. 

Mass polynomials with the symmetry 10, 10, 35, 35 all have factors of (m u — mj). So they only 
appear if we consider the 1 + 1 + 1 case of symmetry breaking. At present we are only considering 
the 2 + 1 case, m u = m c i / m s so we can neglect the 10, 10, 35, 35 representations. 

We found just two singlet tensors in the expansion of 8 <g) 8 <g) 8, so at the symmetric point there 
are only two independent coefficients (usually called F, D or /, d) needed to completely specify all 
the matrix elements between the members of the octet. These give the classic SU (3) inter-relations 
between octet amplitudes. These are generally found to work rather well. We should however be 
able to do better by also including higher terms in the mass expansion. 

There are 8 octets in the expansion of 8 ® 8 ® 8, so if we work to first order 8m q , the SU (3) 
flavour violation, we have 8 new coefficients. There are still many fewer coefficients than there are 
amplitudes, so there are numerous constraints and cross-relations between amplitudes. The singlet 
and octet tensors are given explicitly in Table 2. This table gives the amplitudes for the baryons p, 
A , £ + , S°; the amplitudes for the other baryons can be deduced from isospin symmetry (which 
we are, for now, treating as unbroken). We have used the notation for the matrix element transition 
B -> B' of 

A tMB = (B'\M\B), (2.8) 

where M is the appropriate operator from Table 1 . We illustrate how this table works by reading 
off the mass expansion for the first two amplitudes 

(pIhIp) = a nt]N = V3f-d + (n-s 2 )Smi 
(£ + |tj|I+) =A tllZ = 2d + ( ri +2V3r 3 )5mi 

': (2.9) 

3. 'Fan' Plots 

In the case of hadron masses we found that 'fan' plots were a useful way to display our results, 
see e.g. Fig. 20 of [2]. We plotted the masses of the hadrons in a multiplet against 5m/ = mi — m. 
When 8mi = (the SU(3) symmetric point) all masses are equal, as we increase the symmetry 
breaking the masses fan out, about an 'average' mass which is almost constant. 

We can display matrix elements in a similar way. Some matrix elements, however, will be 
protected from first order flavour symmetry breaking effects. The Ademollo-Gatto theorem, [6], 
for example, states that certain form factors for vector currents will not display 8mj effects. 
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1 a b'mb 


1 

/ ^ 


8 

ri r 2 r 3 r 4 r 5 si s 2 s 3 


NriN 
EtjE 
At]A 
St]S 


\/3 -1 
2 
-2 
— V3 -l 


l -l 
l 2\/3 
l 2 
1 1 


1 NkN 
1 l7rl 
1 SrcS 


1 \/3 

2 
l -V3 


00 -2 20 
00 00 -2 ^3 
00 20 20 


1 ZjtA 
1 AnL 


2 
2 


1 -V3 i 
1 -V3 -i 


\ NKL 
\ NKA 

\ ^ KE 
\ ^ KE 


-V2 \/6 
-V3 -l 
\/3 -l 
y/2 yfl 


y/2 iy/2 y/2 iy/6 
1 i iV3 -a/3 1 — / 
1 -i -/\/3 \/3 -1 -i 
\/2 -/\/2 -\/2 /\/6 


\ AKN 

\ ^ A 
\ SAT 


-y/2 yfl 
-\/3 -l 
V3 -l 
V2 yfl 


y/2 -/\/2 \/2 -/\/6 
1 -i -iy/3 -\/3 1 i 
1 i iy/3 V3 -1 i 
\/2 iV5 -\/2 -i\/6 



Table 2: Coefficients in the mass Taylor expansion of operator amplitudes: 51/(3) singlet and octet. These 
coefficients are sufficient for the linear expansion of hadronic amplitudes. 



3.1 The /-fan 

Using Table 2 we can construct five quantities F{, which all have the same value (2/) at the 
symmetric point, but which can differ once SU (3) is broken. 



1 2 

Ft = —(A^-Agjjg) = 2f- -^s 2 Smi 

?2 = OW+Ag^) = 2f + 4s l 8mi 

^3 = A tnZ = 2f+ (-2fi + V3s 2 )8m, (3.1) 

F 4 = ^=R(A iKZ -A NKZ ) = 2f-2s x 8m l 

1 2 

F 5 = -^(A^s -A RKA ) = 2f+-—(V3s 1 -s 2 )8mi. 

Plotting these quantities gives a 'fan' plot with 5 lines, but only 2 slope parameters (s\, s 2 ), so the 
splittings between these observables are highly constrained. 

A useful 'average F' can be constructed from the diagonal amplitudes 

X F = -(3F 1 +F 2 + 2F 3 )=2f + 0(8mj). (3.2) 
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We expect that a 'fan' plot of F; = Fi/Xp might be less noisy than a plot of F; alone. (Using F; 
rather than F; would also remove renormalisation constants, see eq. (4.10).) In general we shall 
denote quantities with a tilde that have been normalised with an appropriate X. 

3.2 The J-fan 

Similarly, we can construct seven quantities D,, which all have the same value (2d) at the 
symmetric point, but which can differ once SU (3) is broken. 

£>i = -(Ajv^+Agjjs) = 2d — 2r\8mi 

D 2 = A^y. = 2d+ (n + 2V3r^)dmi 
D 3 = -A Ar)A = 2d - (n+2r 2 )8mi 
1 4 
° 4 ~ ^ A ^ N ~ A ^ = 2d ~-^ r 3 Sm i ( 3 - 3 ) 

D 5 = 9U £?rA = 2d+ (r 2 - V3r 3 )8 mi 
1 2 

° 6 ~ 7E^ Ankz +AtK ^ = 2d+ r ^ 5mi 

D 7 = -K(A mA +A Mz ) = 2d-2r 2 5m l . 

Plotting these quantities gives a 'fan' plot with 7 lines, but only 3 slope parameters (r\,r 2 and r$), 
so once again the splittings between these observables are highly constrained. Again it is possible 
to construct an 'average D' similar to 

X D = -(D l +2D 2 +D 4 )=2d + 0(8mf), (3.4) 
in order to produce a less noisy 'fan' plot. 

4. Lattice Calculations 

In order to extract the matrix elements for some operator O, it is necessary to take an appro- 
priate ratio of three and two-point correlation functions, [7, 8] 

c^ B \t,v,p,p>) \ (*(t,pr)c»{x,pr)c*{t-%,p) 

cf(t, P >) V c»( hP )cB(T,p)c?(t-<T,p>)' K - } 

where 



qr» (t,T,p,i/) = tx D r(3g B ,(t;p')O(x;q)@ B (0;p)) 

unpol 



C B {t,p) = tr D r unpol (& B (t;p)& B (0;p)). (4.2) 



This is designed so that any smearing for the source (at time 0) and sink operators (at time t) is 
cancelled in the ratios, [9]; of course smearing improves the overlap with the lowest lying state. 
The baryon operators used are as follows 



(t;p) = £e ip -*e ijk u ia (x,t)(uf(x,t)Cy5d k (x,t)) 
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&La(f,p) = Y, e ip ' X £ijkU ia (x,t)(uf(x,t)CY 5 Sk(x,t)) 

x 

■@Aa(f,p) = Y, e ~ iP '* e ijkSia{x,t)(uf(x,t)CY 5 d k (x,t)) 

X 

&Za(t',p) =^e- ip % k s ia (x,t){sf(x,t)Cy 5 u k {x,t)), (4.3) 

x 

where C is the charge conjugation matrix and i,j,k are colour indices and a is a Dirac index. The 
u and d quarks are treated distinctly, but with degenerate mass. The transferred momentum from 
the initial, B, to final, B' state is given by 

q = p-p'=(i(E B (p)-E B ,(p')),p-p') . (4.4) 

In this study we shall restrict ourselves to zero 3-momentum transfer, i.e. p = p' = when 

q 2 ^ql ax = -(M B -M B ,)\ (4.5) 

such that postive q 2 is spacelike while negative q 2 is a timelike quantity. The energy of the initial 
and final states are now simply the rest masses M B and M B i respectively. We shall also consider 
the time-like component of the vector and axial-vector currents as described in Table 1 by taking 
M = 71° (or 7T + ) and y = y 4 or M = K + and 7 = 7375 respectively. For example we can set O to be 

V 4 = u Y4 d 1 for ^ = Q V 4 = u Y4 s 1 for ^ = 1 ^ 

M = uj 3 y 5 d J A 3 = uj375S J 

Even though in this case p = p' = 0, in general there is still a 4-momentum transfer and we usually 
have a non-forward matrix element. Thus we have 

^^ A B 'MB(l 2 max), *S t ,t — T > °° , (4.7) 

depending on whether we are considering the vector O = V4, F = F unpo \ or axial- vector O = A3, 
r = Fp i three-point function. 

The computed matrix elements are bare (or lattice) quantities and must be renormalised, 

V 4 R =Z V V 4 , A R 3 =Z A A 3 , (4.8) 

where we have denoted the renormalised matrix elements with a superscript R. If is known 

then the renormalisation constant can simply be determined from 

A- A- 

7 _ B'M V B 7 _ B'M A B ( a q\ 

Zy — , z A — . (4.y) 

a b'm v b a b'm a b 

Alternatively by considering ratios in the 'fan' plots the renormalisation constant cancels, for ex- 
ample 

F A = ^ = ^=F A «, i = l...,5, (4.10) 
and similarly for py , £)f and Dj . 
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The renormalised vector / and d coefficients are known (i.e. the matrix elements at the SU (3) 
flavour symmetric point). The vector current, being conserved there, essentially just counts the 
number of quarks. For example using eq. (2.9) we have 

A %vn = -^(2+l-0) = V3^-4, 

= ^=(2 + 0-2) = 24, (4.11) 

giving 

fv = ^=, and 4 = 0. (4.12) 

This result can be used to estimate the renormalisation constant. Either eq. (4.9) can be used at the 
symmetric point or equivalently if we have measured Xp v (7 = 74) then from eq. (3.2) we have up 
to 0(8mf) 

A' ------ , (4.13) 

A F V Ap v Ap v 

For the axial current we can connect our conventions with others in the literature, e.g. [12, 13], via 

f A = ±=F A , and d R A = ^=D A , (4.14) 

at the symmetric point. Similarly to Zy for Z A again either eq. (4.9) can be used, for example for 
neutron j3 -decay 

A NKN = ^8A, (4-15) 

where g A the axial- vector coupling in /3 -decay at the physical point. Equivalently from measuring 
Xf A and Xj) A yields 

y X F A VI y X D A ld R A 

a = = > or Z A = — = — , (4.16) 
Fa a f a *d a *d a 

(provided that either Xp or Xfc is known). Alternatively the ratio where Z A cancels is 

4 = ^L. (4.17) 

d A Xd a 

5. Results 

In a similar manner to previous simulations [1,2] our gauge field configurations have been 
generated with Nf = 2 + 1 flavours of dynamical fermions, using the tree-level Symanzik improved 
gluon action and nonperturbatively 0(a) improved Wilson fermions, [14]. The quark masses are 
chosen by first finding the St/ (3) flavour symmetric point where flavour singlet quantities take 
on their physical values and vary the individual quark masses while keeping the singlet quark 
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Ensemble 






1 


0.12083 


0.12104 


2 


0.12090 


0.12090 


3 


0.12095 


0.12080 


4 


0.12100 


0.12070 


5 


0.12104 


0.12062 



Table 3: Ensembles used in calculations here. 



mass m = (m u + md + m s )/3 = (2m; +m s )/3 constant, as described in section 2. Simulations 
are performed on lattice volumes of 24 3 x 48 at /3 =5.50 corresponding to a lattice spacing of 
a « 0.079, [2]. Our calculations are performed on 5 ensembles chosen by methods as also outlined 
in [2] presently using ~ 400-500 trajectories for off-diagonal matrix elements B' ^ B and ~ 2000 
trajectories for diagonal matrix elements B' = B. In Table 3 we give these (jq, K s ) values used 
here. Ensemble 2 corresponds to the symmetric point. Note that for ensemble 1 we have a universe 
where the / quarks are heavier than the s quark. 

In [5] (see also [2]) the value of the distance away from the symmetric point to the physical 
point was determined to be Strf = —0.01102(3) in lattice units. (The * denotes the physical point; 
see eq. (2.3) for the definition of 5m/.) 



5.1 Xd, Xp 

We first consider Xp, eq. (3.2), and Xd, eq. (3.4) for the axial-vector case. In Fig. 2 we plot 



0.7 r 

0.68 
0.66 



0.64 - 
j 0.62 - 
0.6 - 



X 0.58 - 
0.56 - 
0.54 - 
0.52 - 



1" 



0.5 > 1 1 1 1 1 ' 1 ' 1 

-0.012 -0.01 -0.008 -0.006 -0.004 -0.002 



x, - 



0.002 0.004 0.006 



Figure 2: X f a (filled squares), Xqa (filled circles) against 5m/ each together with a constant fit. 

XpA , Xjy, against 5m/. We expect these quantities to be constant (up to 0(8mf) terms) and within 
error bars this is indeed the case. From the values of the constant fits and using eq. (4.17), we find 
/a/cIa = 1.08(3). As it is well known that axial current matrix elements suffer from large finite 
size effects, we are presently repeating the determination of the matrix elements on larger 32 3 x 64 
lattices. 
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5.2 'Fan' plots 

We now turn to a discussion of 'fan' plots. As noted previously [2], it is better to consider 
ratios, which are less noisy. In Fig. 3 we show F A = F; /Xp A (left panel) and tjf = Df /Xjj a (right 




5 m ! 5m j 



Figure 3: Left panel: F, A = Ff /X Fa (i = 1 , . . . , 5), right panel: Bf = Df /X Da (/ = 1 , . . . , 7). Shown are just 
the numerical results for F A (circles), F A (squares) and &\ (circles), (squares) against dmi together with 
the (normalised) fits from eqs. (3.4), (3.2). 

panel) together with representive numerical results. (Due to our relatively low number of con- 
figurations used in the analysis, error bars overlap if all numerical results are plotted.) Note that 
because we have normalised the data with X f a or Xqa so at the symmetric point the ratios are 1 
exactly. 

We make a simultaneous fit to the data to arrive at a determination for the constants sf, s%, 
and rf, r^, (where sf = sf/2f A , ff = rf/2d ) based upon eqs. (3.2), (3.4). Using these parameters 
it is then possible to reconstruct the equations describing the SU (3) flavour symmetry breaking 
effects on matrix elements up to 0(8m{). It should be noted that we have yet to perform calculations 
for certain correlators, though we are progressing in this direction [10]. Those missing include 
Azk a e wmc h then precludes the use of F A and in our simultaneous fits. The axial-vector matrix 
elements clearly have linear terms in 8m\. 

For the vector case due to the Ademello-Gatto theorem, [6], the linear terms in 8m\ are absent. 
So instead of a 'fan' plot we show in Fig. 4 a selection of matrix elements Aj^ K v A , A A7tV £, and 
^ak v ~- ^ s a check, from Table 2 we see that to leading order A^- v = + d v i.e. it contains only a 
dy term and no f v term. However d v = 0, eq. (4. 12), so we expect to 0(8mj) that A^j vanishes, 
which is clearly seen in Fig. 4. The other decays are also flat. At the symmetric point (or indeed at 
other points) we can estimate Zy from eq. (4.9). As A^ RVA = — a/3/2 = —A^y^ then we expect 
one result in Fig. (4) to be the mirror image of the other. This is the case and choosing Aj^ K v A gives 
upon using eq. (4.9), Zy = 0.87(5). 

6. Octet hyperon semi-leptonic decays 

The theory outlined in previous sections is general; most phenomenological calculations are 
directed towards the semi-leptonic decays B — > B'eV e of various octet hyperons in order to help 
determine \V US \, e.g. [11]. We now briefly indicate how far our programme has reached this goal. 
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2 

1.75 
1.5 
1.25 
1 

0.75 
0.5 
, 0.25 
! 
-0.25 
-0.5 
-0.75 
-1 
-1.25 
-1.5 
-1.75 
-2 



-s— f- 



5 -I I 



i 




-0.012 -0.01 -0.008 -0.006 -0.004 -0.002 0.002 0.004 0.006 

5m, 



Figure 4: The matrix elements Aj^ K v A (red filled circles), A Ak v^ (green filled diamonds) and A-^v? (blue 
filled squares) against 8mi, together with a quadratic fit. 

In the Euclidean metric the general form of the matrix element for semi-leptonic transitions 
B — > B'ev e is 



where 



and 



(B'(p')\Va(q)+A a (q)\B(p)) = a B ip^(q) + ^( q ))u B (p), (6.1) 



ffliq) = JafAq 2 ) + O aB q B f ^ +iq a f ^ , (6.2) 



(Note that in our definition we follow [15], by symmetrising the mass terms appearing in the 
denominator.) The form factors ff (vector), /| (weak magnetism) and /| (induced scalar) cor- 
respond to the vector component, while g" (axial-vector), g^ (weak electricity) and g R (induced 
pseudoscalar) correspond to the axial-vector component of the current. 

Our longer term aim is primarily to determine the CKM matrix element \V US \ from the AS = 1 
semi-leptonic decays, [12, 16] 



Gl 



60k 2 



(M B -M B ,) 5 (l-3S)\V us \ 2 \ff(0)\ 2 1+3 



*?(0) 



/f(0) 



+ ... , (6.4) 



where Gp is the Fermi constant, 8 = (M B — Mg ) / (M B + M B i ) . Hence for a determination of | V us \ 
we require a knowledge of the form factors ff(q 2 ) and g\{q 2 ) at zero 4-momentum transfer, q 2 = 0, 
together with a chiral extrapolation to the physical point. (When at zero momentum transfer the 
form factors ff and g R { are simply the vector, gy, and axial-vector, g\, coupling.) Although this is a 
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complementary determination to the more common kaon semi-leptonic decay determination, it is 
more complicated, not least because it involves an axial form factor. 

However, as a first step, in this work we have restricted our calculations to the specific case 
of zero 3-momentum transfer. Thus from eq. (6.2) for the vector case, we compute the linear 
combination 

Ajfgy B (<7max) = /l (<7max) ~~ _|_ ^ /I (<7max ) = fo (<7max) ■ (6-5) 

Similarly for the axial vector case, from eq. (6.3) we have the combination 

A g'^A B (9max) = £l(4max) ~ ^ + ^ #2 0?max ) = ^o(^max) ■ (6-6) 

(The notation /q, g$ is customary for these form factor combinations.) This is not enough to 
determine ff and g* at q = and at the physical point. (Form factors of matrix elements are 
functions of both q 2 and 5m;.) So in order to disentangle the form factors and to explore the effects 
on ff, gj of symmetry breaking it will be required, in the future, to examine these form factors at 
various values of transferred momenta so that they can be separated as discussed in [15, 17, 18]. 
Phenomenological analyses are given in e.g. [11, 13, 19]. In particular [19] introduces a method 
similar to ours. The Ademello-Gatto theorem actually complicates the determination: we now 
have to find small second order SU (3) flavour symmetry effects, which appear here to be very 
small. 

Note that this determination does not require an explicit determination of the lattice renormal- 
isation constants Zy and Za; as we are interested in deviations from the SU (3) flavour symmetric 
value, it is sufficient to normalise the result either at the symmetric point or (equivalently) with the 
averages Xp v and X^ A , Xp A . 

7. Conclusions 

We have taken the first steps in determining SU (3) symmetry breaking effects for matrix ele- 
ments of all bilinear quark operators for the baryon octet. The strategy of lattice simulations from a 
point on the SU (3) flavour symmetric along the path to the physical point keeping the average quark 
mass constant is ideally suited to investigating these effects. While of intrinsic interest themselves, 
of more phenomenological interest is the determination of form factors relevant to determination of 
the CKM matrix element |V,, V |. This requires more complicated momentum transfer computations, 
and an investigation of 0(8mj) effects, both of which we are now embarking upon. 
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